AD-Al'19  262  ILLINOIS  UNIV  AT  URBANA  DEPT  OF  MECHANICAL  AND  INDUST— ETC  F/G  19/1 

MODELING  OF  DEFLAGRATION-TO-SHOCK-TO-DETONATION  TRANSITION  (DSD--ETC(U) 
JUL  82  H  W  KRIER» . P  B  BULTER*  D  COYNE  AFOSR-81-0145 

UNCLASSIFIED  UlLU-ENG-82-4002  AFOSR-Tr-82-0686  NL 


Annual  Report 

MODELING  OF  DEFLA8RATI0N-T0-SH0CK-T0-DET0NATI0N  TRANSITION 
POROUS  HIGH  ENERGY  SOLID  PROPELLANTS  AND  EXPLOSIVES 


Prepared  by: 

Herman  Krler,  Principal  Investigator 
P.  Barry  Butler  and  Daniel  VI.  Coyne, 
Graduate  Research  Assistants 

Department  of  Mechanical  and  Industrial  Engineering 
144  Mechanical  Engineering  Building 
1206  West  Green  Street 
University  of  Illinois  at  Urbane-Champaign 
Urbana,  IL  61801 


Prepared  for: 

Air  Force  Office  of  Scientific  Research 
Aerospace  Sciences  Directorate 
Bolling  Air  Force  Base,  C  20332 
Progrtm  Manager:  Dr  .  Michael  Strosclo 

For  Research  Performed  under 
Grant  AfOSR  81-0145 
During  Period 

I  dune  1901  through  do  May  1962 


Annual  Report*:  UILU  ENG  82-4002 


MODELING  OF  OEFLAGRATION-TO-SHOCK-TO-DETONATION  TRANSITION  (DSDT) 
IN  POROUS  HIGH  ENERGY  SOLID  PROPELLANTS  AND  EXPLOSIVES 


Prepared  by: 

Herman  Krler,  Principal  Investigator 
P.  Barry  Butler  and  Daniel  W.  Coyne,  Graduate  Research  Assistants 

Department  of  Mechanical  and  Industrial  Engineering 
144  Mechanical  Engineering  Building 
1206  West  Green  Street 
University  of  Illinois  at  Urbana-Champaign 
Urbana,  IL  61801 


Prepared  for: 

Air  Force  Office  of  Scientific  Research 
Aerospace  Sciences  Directorate 
Bolling  Air  Force  Base,  DC  20332 
Program  Manager:  Dr.  Michael  Strosclo 

For  Research  Performed  under 

Grant  AFOSR  81-0145 

During  Period  1  June  1981  through  30  May  1982 


July  1982 


CP IETTTFIC1®SSAH~I  'ATSOI 
A. IK  FOrtCK  0  '  F  ■-  -  :  LE' 

Release;  D1str1bt4fcisitiJo1twifted,78^rt^  wpSR-81,rpi^„rrt  is 

This  tc  ;  ’  V  .  . 

appro*-' '  3 

DistriV.i  t  ;  ‘  " 


matthev  j.  he: 

Chief,  Technical 


Information 


Division 


security  classification  or  this  page  rw*i*n  o«*  e m.r.a; 


REPORT  DOCUMENTATION  PAGE 


REPORT  NUMBER 


AFOSR-TR-  82-06  86 


4.  TITLE  fan4  SuM* Itoj 

MODELING  OF  DEFLAGRATI0N-T0-SH0CK-T0-DET0NATI0N 
TRANSITION  (DSDT)  IN  POROUS  HIGH  ENERGY  SOLID 
PROPELLANTS  AND  EXPLOSIVES 


t.  autho no 

Herman  Krier,  P.  Barry  Bulter  and  Danial  Coyne 


9.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

University  of  Illinois  at  Urbana-Champaign 
Department  of  Mechanical  and  Industrial  Engr. 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1  RECIPIENT'S  CATALOG  NUMBER 


S.  TYPE  or  REPORT  4  PERIOD  COVERED 

ANNUAL  REPORT  -  - 
Julv  1982 


6.  PERFORMING  ORG.  REPORT  NUMBER 

UILU-Eng-82-4002 


«.  CONTRACT  OR  GRANT  NUMBERf*) 


AFOSR  81-0145 


t«.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Office  of  Scientific  Research/NA 

Building  410;  Attn:  M.  Stvoscio 

Bolling  Air  Force  Base,  DC  20332 _ _ 


.  MO* It O RING  AGENCY  HAMS  A  ADDRESS?*'  dilforont  from  Controlling  Oiti co> 


U.  DISTRIBUTION  STATEMENT  (ol  t* U  Rmport) 


10.  PROGRAM  ELEMENT,  PROJECT.  T 
AREA  A  WORK  UNIT  NUMBERS 

Cf(0^  23ol/ 


*2.  REPORT  DATE 

July  1982 


13.  NUMBER  OF  PAGES 


IS.  SECURITY  CLASS,  fof  tM«  raport) 

UNCLASSIFIED 


IS*.  DECLASSIFICATION/ DOWNGRADING 

schedule 


Approved  for  public  release;  distribution  unlimited 


•7.  DISTRIBUTION  STATEMENT  (ol  Uto  mboh-mct  mttorod  in  Mock  20,  It  dfllarcnf  from  J import) 


19.  KEY  WO  ROB  (CauHnuo  on  rover  so  side  ft  nscsmsory  end  Identity  by  block  nimftar) 

Detonation  Physics 
DDT 

j  Porous  explosives 


IMPACT  (Contfmao  on  reverse  H4e  II  mcimiy  and  identity  by  block  nupbtr) 


This  annual  report  represents  the  summary  of  work  done  on  the  modeling  of 
processes  leading  from  deflagration  to  detonation  In  porous  or  granular  high 
energy  propellants.  Particular  attention  Is  paid  to  the  analysis  of  shock  de¬ 
velopment  from  compression  waves  forming  ahead  of  confined  burning  in  the 
original  material.  It  Is  surmized  that  If  the  shock  Is  sufficiently  strong. 
It  will  lead  to  ^shock  to  detonation  ^transition  (SDT).  ^ 


N73 


COITION  OF  t  NOV  49  IS  OBSOLETE 


UNCLASSIFIED 

SECUNITY  CLASSIFICATION  OF  THIS  PAGE  fWfcwi  Dm,  MntmmO 


UNCLASSIFIED 


CCMSiriCATlOW  OF  THIS  PAOCflWnw  Data  Bnfnd) 


MCUniTV 


luring  the  development  of  the  shock  wave,  the  porous  material  may  col¬ 
lapse  into  a  solid  plug  of  void  free  prope^ant  because  the  speed  at  which  the 
wave  propagates  increases  as  the  material  Is  compressed.  The  modeling  effort 
presented  Indicates  how  two-phase  unsteady  combustion  proceses  In  granular  ma¬ 
terial  can  couple  to  the  solid  mechanics  of  shock  formation  and  eventually  to 
a  steady -state  detonation.^ 


Accession  Fo 

~NTIS  C"\JcI 
DTIC  TAB 
Unannounced 
Justif: .atio 

r 

£ 

□ 

n  _  _ 

Distribution/ 

Availability  Codes 

Dist 

A 

Avail 

Spec 

=md/or 

ial 

UNCLASSIFIED 


S1CUHITY  CL  AMiriCATton  OB  PAOCfWh«f*  Ei>CfM/ 


ABSTRACT 


This  annual  report  represents  the  summary  of  work  done  on  the  modeling  of 
processes  leading  from  deflagration  to  detonation  In  porous  or  granular  high 
energy  propellants.  Particular  attention  Is  paid  to  the  analysis  of  shock  de¬ 
velopment  from  compression  waves  forming  ahead  of  confined  burning  In  the 
original  material.  It  Is  surmized  that  if  the  shock  Is  sufficiently  strong. 
It  will  lead  to  ^hock  to  detonation  transition  (SDT). 

During  the  development  of  the  shock  wave,  the  porous  material  may  col¬ 
lapse  Into  a  solid  plug  of  void  free  propellant  because  the  speed  at  which  the 
wave  propagates  Increases  as  the  material  Is  compressed.  The  modeling  effort 
presented  indicates  how  two-phase  unsteady  combustion  proceses  In  granular  ma¬ 
terial  can  couple  to  the  solid  mechanics  of  shock  formation  and  eventually  to 
a  steady-state  detonation. 
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1.  INTRODUCTION 

The  purpose  of  our  work  over  the  past  year  has  been  to  assess  the 
possibility  of  a  detonation  occurring  In  porous  energetic  materials.  In  years 
prior  to  this,  we  were  able  to  show  [1]  that  steady  detonation  characteristics 
(CJ_ pressure,  CJ_ temperature,  detonation  velocity)  were  obtainable  without  the 
stress  wave  analysis  applied  to  the  solid  phase.  This  meant  that  the  rate  of 
flame  propagation  through  the  bed  was  due  entirely  to  the  driving  gas  pressure 
and  not  to  the  shock  initiation  of  the  precursor  wave.  This  was  referred  to 
by  some  as  the  "convective  flame  model." 

In  our  present  work  we  have  developed  a  six-point  outline  which 
summarizes  the  work  in  progress  and  is  a  base  for  future  research.  This 
partitioning  of  the  entire  DSDT  (Deflagration  to  Shock  to  Detonation 
Transition)  process  into  six  parts  has  made  the  more  complex  problem  easier  to 
understand  and  deal  with.  Yhe  logic  followed  in  developing  it  was  as 
follows.  Compression  waves  of  increasing  strength  are  continually  propagating 
ahead  of  the  flame  front  and  eventually  coalesce  into  a  discontinuous  shock 
front.  During  the  development  of  the  shock  wave,  the  porous  material 
collapses  into  a  solid  plug  of  void  free  propellant.  Because  the  speed  at 
which  the  waves  propagates  Increases  as  the  material  is  compressed,  the  shock 
wave  is  at  all  times  increasing  in  strength.  If  this  shock  grows  to  a  large 
enough  strength  within  the  physical  constraint  (length)  of  the  material  being 
shocked,  it  may  cause  detonation  by  one  of  several  methods  to  be  discussed 
later  in  the  report. 
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2.  THE  MODEL 

Shown  In  Fig.  1  Is  a  schematic  representation  of  the  six  part  DSDT. 
Superimposed  on  each  figure  Is  a  line  representing  the  gas  porosity  (gas 
volume/total  volume)  as  a  function  of  bed  location.  A  value  of  <j>  equal  to 
unity  represents  a  zone  of  all  gas  while  $  equal  to  zero  Indicates  a 

homogeneous  solid. 

Part  A  represents  the  Initial  blast  caused  by  the  ignitor  used  to 
Initiate  DSDT  in  the  bed.  In  the  second  portion,  the  convectlvely  burning 
porous  bed  Is  shown  to  be  pressurizing  the  region  behind  the  Ignition  front. 
|This  Is  a  result  of  the  product  gases  being  confined  by  the  walls.  The 

porosity  variation  for  this  time  Indicates  an  all  gas  makeup  In  this 
particular  zone.  Actually,  there  is  a  finite  reaction  zone  where  the 
propellant  Is  decomposing  into  the  product  gases.  Reference  [1]  showed  that 
if  the  reactive  material  was  made  of  large  enough  particles,  this  reaction 
zone  may  be  very  large.  This  results  in  a  slower  pressurization  rate  of  the 
gas  cavity.  The  significance  of  this  on  the  shock  development  will  be  made 
clear  in  Section  D. 

Illustrated  In  the  last  four  parts  of  Fig.  1  is  what  we  refer  to  as  the 

second  possible  branch  of  DSDT,  the  shock  compression  branch  (also  see 

Fig.  2).  This  will  be  expanded  on  In  Sections  C  through  F. 

As  Indicated  by  Fig.  2,  we  believe  that  there  may  be  mechanisms  other 
than  shock  compression  theory.  Because  we  are  dealing  with  a  porous  mixture, 
hot  gases,  as  well  as  stress  waves,  may  be  propagated  ahead  of  the  flame 
front.  This  Is  untrue  when  analyzing  the  case  of  homogeneous  solid  where  only 
stress  waves  can  be  propagation  ahead. 

What  we  refer  to  as  Branch  1  In  Fig.  2  could  occur  If  the  propellant 
grains  are  small  enough  and  their  chemical  composition  energetic  enough.  In 


FIGURE  2.  DSDT  FLOW  CHART 
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this  branch,  the  rate  at  which  the  flame  burns  or  regresses  through  the  open 
propellant  structure  may  have  detonation  characteristics. 

This  rapid  deflagration  is  not  the  same  as  the  steady -state  open  ended 
deflagration  described  In  textbooks  dealing  with  reactive  gas  dynamics.  But 
because  the  ends  of  the  chamber  are  closed  and  the  energetic  propellant 
provides  for  a  rapid  enough  pressure  rise,  the  flame  accelerates  through  the 
first  portion  of  the  bed. 

The  above  paragraphs  have  given  a  brief  description  of  the  DSDT 
phenomenon  as  we  interpret  It.  The  following  sections  will  expand  upon  our 
knowledge  in  each  of  the  six  steps  in  more  detail.  A  short  discussion  on  the 
different  mechanisms  causing  heating  at  a  shock  front  will  then  be 
presented.  The  final  topic  will  be  the  governing  equations  and  how  they  are 
modified  for  each  of  the  six  different  parts. 

2. A  IGNITION  PROCESS 

Our  model  is  a  one-dimensional  analysis  of  a  porous  propellant  bed  of  a 
constant  cross-sectional  area  and  variable  length,  which  attempts  to  model  the 
apparatus  of  some  experiments  [2,3].  A  uniform  initial  porosity  (gas 
volume/total  volume)  of  0.2  to  0.3  Is  used  to  represent  a  packed  sphere 
configuration.  Except  at  the  ignitor  face,  pressures  and  temperatures  are  at 
atmospheric  levels  throughout  the  bed.  For  our  present  work,  the  ignition  is 
assumed  to  occur  at  the  plane  located  at  x  =  0  where  the  gas  and  particle 
temperatures  are  Increased  to  appropriate  levels.  In  our  future  work,  we  plan 
to  Investigate  initial  pressure  disturbances  created  by  the  ignitor  blast  and 
the  effect  they  have  on  the  OSDT  process. 


2.B  CONVECTIVE  FLAME  SPREADING 

Once  the  ignition  process  is  analysed,  the  next  logical  step  is  to  study 
the  convective  flame  spreading  through  the  energetic  material.  For  this  part, 
the  past  work  of  Butler,  Lembeck,  and  Krier  is  used  [1]. 

Once  ignited,  the  decomposition  of  individual  particles  is  modeled  as  a 
pressure-dependent  regressive  burn.  Conservation  of  mass  requires  the  mass  of 
solid  burning  to  be  transfered  into  the  gas  phase  and  conserved  throughout  the 
entire  control  volume.  The  confined  product  gases  cause  the  local  gas 
pressure  to  rapidly  increase  behind  the  ignition  front  and  as  a  res'  of  the 
pressure  gradients  developed,  hot  gases  are  driven  forward  tl  \gh  the 
unreacted  porous  material.  The  voids  In  the  propellant  allow  this  happen 
where  in  a  homogeneous  naterial  it  would  act  to  only  further  c  the 

product  gases. 

Convective  heat  transfer  from  the  product  gases  to  the  unreacted 
propellant  causes  ignition  of  the  energized  propellant;  the  process  repeats 
itself  as  the  gas  pressure  continues  to  increase.  It  is  this  high  pressure 
gas  zone  which  acts  as  a  driving  force  for  the  stress  waves  being  propagated 
into  the  solid  ahead. 

The  effect  of  the  compression  wave  is  shown  in  Fig.  1C.  At  this  time  the 
leading  compression  wave  has  traveled  Into  the  porous  solid  a  distance,  ax.  In 
the  region  between  the  leading  compression  wave  and  the  flame  front  the 
porosity  decreases  becauses  of  the  increasing  strength  of  the  compression 
waves.  It  will  be  shown  In  the  next  section  that  the  porosity  is  directly 
related  to  the  local  stress  by  a  pore  collapse  model  or  a  pressure-solids 
fraction  (P  -  a  where  a  =  [1  -  $]"*)  law* 

The  analysis  of  a  nonhomogeneous  flow  problem  Involves  the  Introduction 
of  a  new  unknown  Into  the  problem,  the  porosity  (or  mass  fraction).  The 
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porosity  of  the  mixture  can  have  any  value  between  unity  (all  gas)  and  zero 
(all  solid).  To  form  closure  to  the  system,  an  additional  independent 
equation  Is  necessay.  A  possible  equation  to  fulfill  this  requirement  would 
be  a  porosity-pressure  relation,  one  particular  equation  was  the  pore  closure 
model  of  Carroll  and  Holt  [4],  The  next  section  will  outline  this  relation. 

2.C  PORE  COLLAPSE 

The  compression  of  the  porous  matrix  was  modeled  by  Carroll  and  Holt  as  a 
hollow  sphere  where  the  inner  and  outer  radii  are  such  that  the  pore  diameter 
and  the  overall  void  fraction  are  those  of  the  porous  mixture  [4].  The  use  of 
such  a  model  to  represent  a  granulated  bed  of  propellant  spheres  (the  inverse 
problem)  has  not  been  clearly  justified  and  the  possible  errors  associated 
with  it  should  be  considered. 

Pore  collapse  occurs  in  three  distinct  phases.  Phase  1  is  the  elastic 
region.  Phase  2  is  the  elastic-plastic  region,  and  Phase  3  is  the  plastic 
region.  Phase  1  applies  until  the  pressure  has  become  large  enough  to  cause 
plastic  deformation  at  the  inner  radius.  The  change  In  porosity  is  very  small 
during  this  stage.  The  second  phase  or  the  elastic-plastic  region  occurs  as 
plastic  deformation  begins  at  the  inner  radius  and  propagates  outward  to  the 
outer  radius.  As  in  the  elastic  region,  the  porosity  change  in  the  elastic- 
plastic  region  is  very  small.  It  is  during  Phase  3,  where  plastic  deformation 
occurs  throughout  the  solid  material,  that  significant  porosity  changes 
occur. 

The  P  -  a  (pressure-a)  relations  were  derived  by  Carroll  and  Holt  by 
solving  the  momentum  equation  with  an  imposed  time-dependent  hydrostatic 
pressure  P(t)  at  the  outer  radius  and  a  zero  pressure  condition  at  the  inner 
radius  [4],  A  volume  preserving  deformation  assumption  is  also  used.  This 
requires  the  solid  to  be  incompressible  during  the  pore  collapse. 
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Modifications  were  made  by  Kooker  and  Anderson  [5]  to  include  a  porosity 
dependent  yield  stress  (i.e. ,  strain  hardening)  and  a  rate  dependent 
resistance  to  pore  collapse  due  to  material  viscosity.  The  equations  for  each 
of  the  three  phases  for  this  dynamic  model  can  be  found  in  Ref.  [5]. 

When  the  spherical  pore  experiences  a  rate  of  pressurization  less 
than  10  kbar/usec,  a  much  simpler  static  pore  collapse  model  can  be  used. 
This  was  also  presented  by  Carroll  and  Holt  [4].  For  this  model,  the  inertia 
and  viscous  terms  in  the  governing  equations  are  removed  and  the  remaining 
equation  is, 

(aQ  -  1  )  1  Y 

p">- •  <’> 

The  left  side  of  Eq.  (1)  represents  the  pressure  applied  externally  to 
the  sphere.  The  term  Pg<^ (a )  is  the  equilibrium  stress  field  for  the  material 
and  Is  different  for  each  of  the  three  phases  of  compaction.  The  last  term  is 
an  internal  resistance  due  to  isentroplc  compression  of  the  trapped  gases. 
Gamma  (y)  is  the  specific  heat  ratio  of  the  trapped  gas. 

A  further  assumption  is  made  that  no  change  in  porosity  occurs  during  the 
elastic  phase  or  the  elastic-plastic  phase.  This  was  shown  to  be  true  when 
the  ratio  of  yield  strength  to  shear  modulus  is  small  and  the  Initial 
alpha,  aQ ,  is  not  close  to  that  of  all  solid  (a  =  1,0). 

With  this  assumption,  the  P-a  relation  reduces  to  a  single  equation  for 
the  plastic  phase: 

P(t)  =  2/3Y)„  *>  ♦  Pv0  [-^jj  T  (2) 

where  the  transition  pressure  is  approximated  to  be: 
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and  Y  =  yield  stress,  a  =  total  volume'solid  volume  =  (1  - 

In  summary,  a  *  aQ  for  (0  <  P  <  Pcr1t),  a  =  f(P(t)),  by  Eq.  (2)  for 
(P  >  Pcrit)  with  the  following  assumptions, 

(1)  A  porous  bed  can  be  modeled  by  hollow  spheres; 

(2)  The  solid  material  is  incompressible  during  the  pore  collapse; 

(3)  Very  little  porosity  changes  occur  during  the  elastic  and  elastic- 
plastic  phases;  and 

(4)  Dynamic  effects  can  be  neglected. 

Figure  3  shows  the  static  model  for  the  pore-collapse  with  initial  alphas 
of  1.2,  1.4,  1.6,  and  1.8.  The  elastic  and  elastic-plastic  transition 

pressures  are  shown  with  the  dashed  lines. 

The  pressure  required  for  complete  collapse  can  be  seen  from  Fig.  3  to  be 
relatively  small  compared  to  the  pressure  generated  in  the  burning  zone;  thus, 
it  is  anticipated  that  complete  pore  collapse  will  occur  when  the  compression 
waves  reach  a  moderate  strength  (1  to  2  kbar).  This  leads  to  our  fourth  step. 
Plug  Formation. 


2.D  PLUG  FORMATION 

As  stated  in  Section  2.C,  the  particle  stress  to  fully  compact  the 

explosive,  HMX,  is  of  the  order  of  several  thousand  atmospheres.  A  gas 

pressure  of  this  magnitude  will  be  attained  in  a  very  short  time  In  the 

confined  zone  behind  the  flame  front.  This  critical  stress  is  transferred  to 
the  solid  (constant  pressure  boundary  condition)  and  transmitted  through  the 
bed  at  the  local  sound  velocity  plus  particle  velocity.  The  bed  then 


£ 


11 


collapses  as  the  wave  passes  over,  starting  at  the  flame  front  and  progressing 
forward  (away  from  the  products  zone),  forming  the  solid  plug.  This  was  shown 
in  Fig.  Id.  Since  the  high  pressure  gas  zone  can  no  longer  drive  hot  gases 
through  the  permeable  solid,  it  pressurizes  at  even  a  higher  rate. 

For  our  present  analysis,  the  solid  material  is  assumed  incompressible 
during  the  pore  collapse  process.  In  addition,  we  do  not  treat  any  type  of 
frictional  heating  associated  with  the  pore-collapse.  Shear  heating  may  be  an 
important  aspect  and  we  Intend  to  investigate  It  in  the  year  ahead.  When  we 
do,  it  will  probably  yield  a  rate  of  pressurization  dependency. 


2.E  SHOCK  FORMATION 

At  this  point,  the  stress  waves  which  formed  the  solid  plug  have 
coalesced  Into  a  shock  front  traveling  through  the  porous  medium.  This  is 
analagous  to  the  piston  problem  where  stress  waves  sent  out  by  a  piston 
traveling  In  a  closed  tube  eventually  overtake  the  lead  wave  and  develop  a 
discontinuity  or  shock  front. 

The  shock  In  the  porous  bed  forms  a  boundary  between  two  distinct 
zones.  In  front  of  the  shock,  the  material  is  undisturbed  from  its  Initial 
state.  Behind  the  shock,  the  pores  have  collapsed  and  the  solid  matrix  Is 
compressed  to  some  final  state  (Pb>vb).  The  pore  collapse  and  solid 
compression  occur  within  the  thickness  of  the  shock  wave  and  are,  therefore, 
assumed  to  occur  Instantaneously.  This  process  is  discussed  below. 

In  addition  to  the  two  zones  mentioned  above,  there  is  still  the  third 
(product  gas)  zone  Immediately  behind  the  plugged  zone.  Since  the  product 
gases  are  now  completely  confined  (due  to  the  wall  and  solid  plug)  they 
continually  transmit  stronger  stress  waves  over  time,  which  act  to  reinforce 
the  already  developed  shock  wave  traveling  ahead  of  it.  The  reason  for  the 


12 


reinforcement  of  the  shock  wave,  as  in  the  piston  problem,  is  the  property  of 
the  medium  to  Increase  the  local  sound  velocity  as  it  is  compressed.  Because 
of  this,  each  successive  stress  wave  travels  at  a  faster  speed  than  the  one 
prior  to  it,  and  eventually,  they  reach  the  shock  and  strengthen  it. 


2.F.1  SHOCK  INITIATION  OF  POROUS  MATERIALS 

Past  experiments  have  shown  that  most  solids  obey  a  linear  shock 
velocity-particle  velocity  (U$,u)  relation  for  the  shock  strengths  of  interest 
to  us.  That  is,  the  relation. 


Us  =  A  +  Bu 


(4) 


is  a  fairly  good  description  of  experimental  data  available  for  most 
materials.  Here,  A  and  B  are  constants  dependent  on  the  material  being 
shocked.  In  the  limit  of  an  infinitesimal  pressure  disturbance  (sound  wave), 
the  particle  velocity  (u)  approaches  zero  and  the  constant  A  is  then  found  to 
be  the  adiabatic  sound  velocity,  cQ.  The  constant  B  has  no  physical 
significance  other  than  the  slope  In  the  Us  -  u  plane. 

The  linear  form  of  Eq.  (4)  makes  it  convenient  for  developing  other 
important  relations.  Using  it  in  conjunction  with  the  shock  Hugoniot  jump 
conditions,  one  obtains  an  equation  relating  Hugoniot  pressure,  PH,  to 
specific  volume,  v,  namely. 


PH  = 


A2R 


v0(l  -  AR)‘ 


(5) 


where  R  *  1  -  v/vQ.  In  Eq.  (5),  the  subscript  H  is  In  reference  to  the  end 
states  along  the  Hugoniot  and  the  subscript  o  represents  the  initial 
conditions,  also  on  the  Hugoniot  line. 
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Reference  [6]  is  a  compilation  of  shock  data  for  numerous  Inert  and 
reactive  solids.  In  this  reference,  the  U^-u  data  for  each  material  Is  fit 
with  a  straight  line,  l.e.  Eq.  (4).  It  should  be  noted  that  while  the  curve 
fits  for  the  Inert  solids  and  liquids  Involve  many  data  points,  some  of  the 
highly  reactive  explosives  (HMX.ROX)  have  only  two  or  three  scattered  points 
through  which  a  straight  line  Is  fit. 

Figure  4  shows  the  effect  that  a  20  percent  change  in  the  slope  of  a  Us-u 
line  has  on  the  P-v  and  T-v  relations,  the  method  for  determining 
T  =  T(v)  will  be  discussed  shortly. 

Equation  (5)  represents  only  the  end  states  obtained  through  a  shock 
process.  In  other  words,  only  those  states  with  the  pair  (P,v)  on  the 
Hugonlot.  In  order  to  represent  all  states  possible  rather  than  just  the  jump 
conditions,  a  more  general  Mie-Gruneissen  state  equation  is  used,  l.e.. 


P(v,E) 


,  a2  o  -k  <y 

vo  (1  -  BR)Z 


(6) 


One  can  see  by  substitution  that  If  the  general  state  equation  (Eq.  (6)) 
Is  evaluated  along  the  Hugonlot  (e  =  eH  =  0.5  *  PH  *  (vq  -  v)),  then  Eq.  (6) 

becomes  the  same  expression  as  Eq.  (5). 

In  Eq.  (6),  r  Is  the  Grunelssen  coefficient  defined  as  r(v)  =  (3P/3e)y. 
Using  r/v  =  constant  appears  to  be  a  reasonable  assumption  for  shock 

compressions  valid  up  to  vQ/v  =  1.50.  This  will  be  important  in  determining 
the  temperature  rise  due  to  shocking  a  material. 

The  process  of  shock  Initiating  a  homogeneous  solid  will  be  described 
with  the  aid  of  Fig.  5.  Here,  the  solid  material  Is  Initially  at  rest, 
represented  as  state  o.  As  the  shock  fully  traverses,  a  point  on  the 


material.  It  is  compressed  to  state  8.  The  path  process  between  states  o  and 


.HUGONI 


V8  F  v0  vc 


SPECIFIC  VOLUME 

FIGURE  5.  Hugoniots  for  unreacted CX*0)  and 
fully  reacted (X-l)  explosives 
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B  is  not  necessarily  along  the  Hugonoit.  The  Hugonlot  is  only  a  locus  of  end 
points  and  not  a  path. 

If  the  material  was  an  inert  solid  or  liquid  and  the  shock  wave  was  only 
a  short  pulse,  followed  by  an  expansion  wave,  the  shocked  inert  material 
(state  B)  would  expand  back  to  the  original  pressure,  PQ,  along  the  isentrope 
shown  and  have  the  end  state  (c).  However,  when  a  reactive  material  is 
shocked  to  state  B,  the  energy  increase  associated  with  the  shock  wave  will 
cause  some  degree  of  chemical  decomposition.  If  the  amount  of  energy  Is 
sufficient  to  support  the  full  decomposition,  the  products  will  reach 
equilibrium  at  state  CJ,  the  Chapman-Jouget  point.  Here,  A,  the  degree  of 
reaction  Is  unity  representing  complete  reaction.  The  process  of  increasing 
the  energy  of  a  homogeneous  solid  through  the  shock  wave  will  be  discussed  in 
the  next  section,  followed  by  a  description  for  treating  porous  materials. 

2.F.2  TEMPERATURE  INCREASE  DUE  TO  SHOCK  WAVES 

The  method  used  for  determining  temperature  increases  due  to  shock  waves 
Is  that  of  Zel'dovich  and  Raizer  [7]  and  involves  the  partitioning  of  energy 
and  pressure  into  three  separate  components. 

The  first  component  of  pressure  is  due  to  the  interatomic  potential  and 
Is  referred  to  as  the  elastic  or  cold  term.  At  zero  degrees  Kelvin  and  one 
atmosphere  pressure,  the  atoms  neither  repel  or  attract  one  another  and  at 
this  point  the  elastic  pressure  Is  zero.  The  volume  occupied  by  the  solid  at 
this  point,  VQC,  Is  sllgthly  less  than  the  volume  occupied  by  the  same  mass  of 
material  at  room  temperature.  Because  the  solid  thermally  expands  when  heated 
from  0°K  to  300°K,  the  solid  actually  has  a  negative  pressure  (attractive 
interatomic  force)  at  ambient  conditions.  As  the  material  Is  compressed, 
either  Isentroplcally  or  irreversibly  In  a  shock-like  process,  the  atoms  repel 
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one  another  and  therefore  contribute  to  the  elastic  pressure  in  a  positive 
sense.  The  elastic  component  of  the  total  energy  is  related  to  the  elastic 
pressure  through  the  work -energy  relation, 

Pc  dv  =  -dec  (7) 

where  the  subscript  c  Is  In  reference  to  the  elastic  or  cold  component  of  the 
total  Internal  energy,  e,  or  total  pressure,  P.  The  name  cold  energy  or 
pressure  has  significance  in  that  the  elastic  pressure-volume  relation  can  be 
shown  to  be  the  0°K  Isotherm  and  also  the  S  =  0  isentrope. 

The  second  and  third  components  of  the  energy  are  due  to  the  atomic  and 
electronic  motions,  respectively.  The  electronic  terms,  efi  and  Pg,  are  only 
of  significance  above  approximately  10,000°K  and  hence  outside  the  regime  of 
the  processes  encountered  here. 

The  thermal  component  of  energy  represents  the  harmonic  oscillations  of 
the  atoms  and  Is  expressed  as 

eT  "  Cv  (T  *  V  +  V 

To  determine  the  thermal  pressure  from  the  thermal  energy,  eT,  the  first  law 
of  thermodynamics  and  a  Maxwell  relation  are  used  to  obtain 

$7>T  ’  T  <lf)»  -  "•  <»> 

This  1i  a  general  equation  applicable  to  both  components,  elastic  and 
thermal.  Since  the  cold  pressure  Is  only  a  function  of  volume  (or  more 
Specifically  Interatomic  distance,  then  Pc  -  Pc(v)  and  Eq.  (9)  gives 
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the  same  expression  as  given  prior  to  this,  Eq.  (7),  In  the  work-energy 
analysis.  For  the  thermal  component. 


apT 

pt  5  V  <n?v 


00) 


Equation  (10)  can  be  expressed  in  terms  of  the  thermodynamic 

•• 

derivative  r(v)  (the  Gruneissen  coefficient)  as 


ap 

where,  as  stated  earlier,  r(v)  »  v(  -rr)„.  The  Gruneissen  coefficient,  as 

<70  V 

expressed  here,  is  a  measure  of  the  change  in  thermal  pressure  to  thermal 
energy. 

Equations  (7),  (9)  and  (11)  are  solved  in  addition  to  the  algebraic 
component  relations. 


e  -  ec  +  eT  (12) 

P  -  Pc  ♦  PT  (13) 

and  the  Hugonlot  jump  conditions, 

«  •  •«  ■  7  pH<,o  -  *> 


< 


(14) 
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to  generate  the  cold  compression  curve  shown  in  Fig.  5.  For  this  particular 
case,  the  ratio  r/v  is  assumed  to  be  constant.  Figure  5  shows  graphically  how 
each  component  of  energy  Is  represented.  The  area  ABDOA  is  the  total  increase 
in  internal  energy  due  to  shock  pressure,  PB.  Of  the  energy  deposited  by  the 
shock,  the  area  under  the  curve  FED  goes  into  the  potential  or  elastic 
energy.  The  remaining  area  ABEFOA  is  the  contribution  to  the  thermal 
component  and  therefore  determines  the  temperature  rise  associated  with  the 
shock. 


3.  POROUS  MATERIALS 


The  analysis  for  determining  the  temperature  rise  in  a  shocked  porous 
material  Is  somewhat  similar  to  that  discussed  In  Section  2.F.2  except  the 
Initial  volume  Is  now  V00,  larger  than  the  homogeneous  initial  volume,  VQ. 
The  additional  Initial  volume  Is  occupied  by  the  gases  internal  to  the  pores, 
usually  assumed  to  be  air  or  more  simply  a  vacuum. 

Since  the  shock  energy  is  essentially  the  product  of  one-half  the  change 
In  volume  of  the  material  times  the  shock  pressure,  a  porous  mateial  shocked 
to  the  same  final  volume  as  a  homogeneous  material  will  have  more  energy 
deposited  to  it  than  the  homogeneous  solid.  The  elastic  curve  is  the  same  for 
both  cases,  so  the  material  which  was  initially  porous  will  be  shocked  to  a 
greater  thermal  energy,  and  thus  higher  temperature  than  the  homogeneous 
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4.  COMPRESSIVE  WAVES 

In  investigating  the  possibility  of  shock  formation  in  a  porous  bed  due 
to  coalescence  of  compression  waves,  we  have  extended  Macek's  concepts  and 
application  of  the  method  of  characteristics  [9].  The  requirement  of  a  high 
pressurization  rate  or  the  need  to  justify  a  proposed  method  to  increase 
surface  area  Is  not  needed  when  considering  porous  propellants,  as  was 
required  to  explain  the  large  exponential  pressure  rise  (used  by  Macek)  In 
order  for  shock  formation  to  occur  In  a  homogeneous  solid.  Large  surface  area 
Is  present  due  to  the  nature  of  a  porous  material  and  the  high  pressurization 
rates  required  to  form  shock  waves  In  a  solid  are  not  needed  to  form  a  shock 
wave  in  a  porous  material. 

The  first  phase  of  this  work  Is  a  mathematical  model  of  a  porous  bed 
subjected  to  an  Imposed  pressure  gradient  at  one  end.  The  method  of 
characteristics  Is  used  (as  Macek  utilized)  to  determine  if  shock  formation 
from  the  coalescence  of  compression  waves  Is  possible  in  a  porous  mixture.  If 
shock  formation  is  found  to  be  possible,  the  rate  of  pressurization  required 
to  cause  shock  formation  within  reasonable  distances  is  determined. 
Figures  6a  and  6b  show  a  sample  case  run  for  aluminum.  A  second  phase  of  this 

study  will  then  be  to  assess  whether  such  a  pressure  rise  Is  possible  within 

the  burning  zone. 

In  addition  to  employing  the  method  of  characteristics,  work  Is  in 
progress  on  a  numerical  solution  of  the  full  hydrodynamic  equations  governing 
the  stress  wave  build-up.  The  system  of  equations  governing  the  flow  is 

solved  by  utilizing  a  Lagrangian  finite  difference  technique  with  the 

necessary  artificial  viscosity  to  handle  the  severe  gradients  associated  with 
the  shock  development. 
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Figure  6a.  Characteristics  diagram  demonstrating  the  shock  formation  in  a 

solid  (Aluminum,  pQ  =  1.6  g/cc,  CQ  =  2.56  mm/usec)  due  to  compression 

waves  generated  by  an  exponential  pressure  rise  (P  *  .OBe*1*). 


X-LOCATION  (MILLIMETERS) 


Figure  6b.  Development  of  shock  wave  (calculated  from  the  characteristics). 
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5.  HOT  SPOTS  AND  FRICTIONAL  HEATING 

Section  2.F  showed  by  using  classical  shock  hydrodynamics  how  the 
temperature  rise  due  to  a  shock  was  calculated  for  a  homogeneous  and  non- 
hoinogeneous  material.  Here,  the  energy  associated  with  the  shock  wave  was 
balanced  by  the  interatomic  potential  and  thermal  energy  of  the  solid.  The 
porous  material  was  shown  to  have  a  greater  increase  in  temperature  for  a 
given  shock  strength. 

Since  experiments  have  shown  that  explosives  with  internal  voids  shock 
initiate  at  lower  Impact  velocities  (lower  shock  pressures)  than  homogeneous 
explosives,  the  Idea  of  localized  heating  or  "hot  spots"  is  a  topic  of  current 
Interest.  This  means  that  there  may  be  additional  mechanisms  to  explain  the 
initiation  of  porous  explosives  beyond  the  rise  in  bulk  temperature  due  to  its 
porous  nature. 

Hayes  [8]  presents  a  model  for  preferentially  depositing  a  fraction  of 
the  shock  energy  to  a  certain  portion  of  the  material  (proportional  to  the 
initial  void  fraction)  and  depositing  the  remaining  energy  to  the  bulk 
material.  The  bulk  material  Is  assumed  to  be  compressed  isentropical ly  while 
the  hot  spots  are  compressed  irreversibly.  This  leads  to  two  distinct 
temperatures  in  the  material. 

Frey  [10]  focuses  on  the  heating  of  porous  materials  due  to  friction  of 
particular  Interest  to  us  in  his  extension  of  the  work  by  Carroll  and  Holt  to 
determine  the  temperature  rise  at  the  interior  of  the  pores  at  different 
pressure  rise  times.  This  analysis  is  important  in  determining  whether  or  not 
a  reactive  material  will  ignite  upon  pore  closure. 
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6.  CONCLUDING  COMMENTS 

This  brief  report  has  provided  Information  on  the  Inter-connected,  fluid 
and  solid  mechanical  processes  that  must  be  described  if  one  is  to  understand 
the  mechanism  for  deflagration  to  detonation  transition  (DOT)  in  porous 
propellant  explosives.  It  is  clear  that  there  remains  considerable  work  to  be 
done,  not  the  least  of  which  is  to  solve  the  unsteady  shock  reactive  flow 
conservation  equations.  An  extensive  "hydro-code'1  is  being  developed  to 
numerically  solve  the  battery  of  equations  which  contains  the  system. 

Whenever  possible,  information  from  the  latest  reported  data  is  utilized 
to  formulate  the  many  constitutive  relations  required  for  closure.  We  expect 
to  present  useful  predictions  for  DDT  (and  DSDT)  in  porous  reactive  solids 
based  on  work  continuing  on  our  continued  AFOSR  supported  work. 
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